Numerical Relativity meets Data Analysis: 
Spinning Binary Black Hole Case 



> 



OO ■ Deirdre Shoemaker, Birjoo Vaishnav, Ian Hinder and Prank 

^ ■ Herrmann 

' Center for Gravitational Wave Physics, Penn State, University Park, PA 16802 

D 

I ■ Abstract. We present a study of the gravitational waveforms from a series of 

spinning, equal-mass black hole binaries focusing on the harmonic content of the waves 
' and the contribution of the individual harmonics to the signal-to-noise ratio. The 

gravitational waves were produced from two series of evolutions with black holes of 
^ ' initial spins equal in magnitude and anti-aligned with each other. In one series the 

magnitude of the spin is varied; while in the second, the initial angle between the black- 
$-H I hole spins and the orbital angular momentum varies. We also conduct a preliminary 

investigation into using these waveforms as templates for detecting spinning binary 
black holes. Since these runs are relativity short, containing about two to three orbits, 
merger and ringdown, we limit our study to systems of total mass > SOMq. This choice 
ensures that our waveforms are present in the ground-based detector band without 
Ps| i needing addition gravitational wave cycles. We find that while the mode contribution 

' to the signal-to-noise ratio varies with the initial angle, the total mass of the system 

. caused greater variations in the match. 
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1. Introduction 

Gravitational waves produced during the coalescence of compact objects such as black 
holes, are one of the most promising sources for detection by interferometric detectors 
such as LIGO, VIRGO, GEO and TAMA [IJ. For these ground based-detectors, where 
any gravitational-wave signals may be deeply buried in detector noise, the matched 
filtering technique [2] is the detection strategy of choice. Matched filtering is optimal 
when accurate representations of the expected signal are used. For low mass compact 
object binaries that means templates built from well-known analytic methods such 
as the post-Newtonian (PN) approximation [3J. When the total mass is larger than 
approximately SOMq [3], the merger of the binary black hole (BBH) system will not 
only be present in the sensitivity band of ground-based detectors but also generate the 
strongest signal. Quantifying the last orbits and merger of a BBH system has long been 
the purview of numerical relativity. 

A new direction in the efforts to generate accurate templates is the inclusion of 
waveforms produced by numerical relativity. The waveforms from numerical relativity 
may be used in data analysis methods involving searches for inspiraling sources 
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in several ways. For instance, they can be used in creating hybrid templates for 
detection, validating and extending the mass range of the current search strategy and 
may be implemented as templates themselves especially for the higher mass range 

[HlEllTlElElIlolin]. 

In a previous paper [9] (paper I), we studied waveforms from a series of numerical 
evolutions of equal-mass, spinning BBH mergers, which we label the A-series. In the A- 
series, the black holes had spins that were equal in magnitude with one spin aligned and 
the other anti-aligned with the orbital angular momentum, J^rb- We found that using 
only the dominant harmonic of the radiation, a common practice in numerical relativity 
and data analysis, caused a complete degeneracy of the this spin space with respect 
to detection. In other words, a template of a non-spinning BBH waveform faithfully 
matched the signal of any anti-aligned spinning waveform when only the dominant modes 
were used. However, retaining non-dominant harmonics of non-negligible amplitude 
broke the degeneracy to some degree, most notably for spins J/M^ > 0.6. In this 
paper, we further our previous study in two ways. First we include a new series of BBH 
mergers, the B-series, in which the initial black-hole spins are still equal in magnitude 
and anti-aligned with each other but the initial angle the black-hole spins make with the 
orbital angular momentum, is allowed to vary. Second, for both the A and B series 
of data we conduct a new study of the contributions of the individual modes to the 
signal-to-noise ratio. In addition, for the B-series we also calculate the dependence of 
the minimax matches between the dominant mode and the full waveform on the initial 
spin orientation. 

2. The Waveforms 

The binary black hole coalescence problem can be divided into three phases called the 
inspiral, merger and ringdown. Due to the early accessibility of PN and other analytic 
based approaches, most of the work in setting up detection schemes has been done 
with these analytic approaches. These methods are well suited for the inspiral phase of 
the coalescence, only breaking down at some yet-to-be-determined point within a few 
orbits of the merger. Fortunately, the signal resulting from the inspiral is in the ground- 
based frequency band for systems of total mass less than approximately SOMq. As the 
mass increases, the inspiral lowers in frequency, and the detectable signal contains more 
merger and ringdown. Now that numerical relativity is producing the waveforms for the 
final orbits and the merger phase of the coalescence IT3l E] we can investigate the 
inclusion of the merger regime in detection strategies. 

We study two sets of waveforms both the result of evolutions conducted by the PSU 
numerical relativity group. The A-series was published in [15] with initial black-hole 
spins covering the set a = 0.0, 0.2, 0.4, 0.6, 0.8. The B-series was published in [16] and 
generalizes the A-series with variation of the initial angle that one of the anti-aligned 
spins makes with the axis, i), at a fixed magnitude of spins, a = 0.6. When = 0, we 
recover the a=0.6 waveform of A-series. When ■(9 = 7r/2, the spin-directions lie in the 



Spinning Binary Black Holes 



3 



plane of the orbit and are in the "superkick" [T7j configuration in which the maximum 
gravitational recoil from the BBH mergers has been found. 

Since these waveforms were originally produced to study the gravitational recoil 
imparted to the final black hole after an asymmetric collision, only two to three orbits 
were evolved (the merger phase dominates the recoil). The number of orbits is set by 
the initial orbital frequency for a given total mass. In order to place the numerical 
waveforms firmly in the frequency band of the detector, we use the initial LIGO noise 
curve [IE] and we only investigate masses larger than SOM© when calculating matches 
between waveforms. The total mass sets the frequency at which the signal enters the 
band. For example, the cutoff frequency for a binary system of SOM© is 0.02/M or 
approximately 80Hz and about 40Hz for IOOMq. 

The waveforms were extracted from the numerical evolution of the spacetime in 
terms of the Newman-Penrose scalar, /^{t,x,y, z), which is expanded into angular 
modes via _2l£m(6', 0), the spin-weighted s = —2 spherical harmonics, by extraction 
on a sphere. The dominant mode for the quasi-circular orbits is the quadrupole mode 
= m = 2). The angles 9 and (f) correspond to the inclination and azimuthal angles 
between the source and detector in the source frame. When ^ = 0, the observer is 
directly above the orbital plane of the binary and sees primarily the i = m = 2 mode. 
As 6 increases, the waveforms are a mixture of modes. We truncate the infinite series of 
modes at £ < 4 because we do not extract all the modes from the simulations and modes 
of £ > 4 were zero within our numerical error. As more complicated configurations are 
evolved, more modes will need to be accurately extracted from the codes. 

3. Faithfulness 

The multipolar analysis of BBH waveforms produced by numerical relativity has been 
pursued for both unequal mass and spinning BBH configurations [191 1201 [21|. In paper 
I, we found that using only the dominant mode in comparing waveforms, tantamount 
to choosing an inclination angle with the detector in the source frame of 6* = 0, resulted 
in a degeneracy of the A-series parameter space. In this paper, we focus our attention 
on how different initial configurations, in this case a and i), result in different mode 
contributions to the signal-to-noise ratio. To build intuition about what parameters 
might be important to the template space of black-hole mergers, we further calculate 
the overlap between pairs of our waveforms. We perform a preliminary analysis on 
how faithful i = m = 2 waveforms of various parameters would be in matching with 
waveforms at random inclination angle. We keep the total mass for each template fixed 
and vary the inclination angle of the detector, 6, the spin a, and initial angle when 
appropriate. 

The minimax match is given by [22l [23] 



M = max min 

to <3> 
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where 



('.IM = 4Re/:7»^^./. (2) 



The Fourier transform of the strain, h{t), is given by h+{f) = jF(Re(^'4)) (/)/(— 47r^/^) 
where \l/4(t) = ^{h+{t) — ihx{t)) , and JF is a Fast Fourier Transform. The signal-to- 
noise ratio, p, is given by 



[ "'/min Shif) J 

The variable Sh{f) denotes the noise spectrum for which we use the initial LIGO noise 
curve. The domain [/min, /max] is determined by the detector bandwidth and the masses 
of our signal. The masses are set such that the overlaps will not change significantly 
if we were to add the inspiral portion of the signal because the A and B series of 
waveforms have orbital frequencies that increase almost monotonically. Owing to this, 
the gravitational wave frequency also increases monotonically with time implying that 
extending the signal back in time will not change the spectrum in the merger band. 
When precession is significant this will no longer be true and the inspiral will likely 
contribute to the signal at higher frequencies. 

A-series: In Fig. [T]we plot the match versus spin at different inclination angles, 6, for a 
given total mass of lOOM©. This plot first appeared in paper I and is included here for 
reference. The figure shows that as the spin increases, the match between a waveform of 
^ = and one of non-zero 6 decreases. This indicates that the non-dominant modes are 
important both for distinguishing between different spinning waveforms and in making 
a detection. This is most notable for the a = 0.8 case. 
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Figure 1. The minimax match versus a for three values of 9, 7r/4, 7r/3, 7r/2 for the 
case of total mass IOOM0. The value of = is not included in the plot, since this is 
the waveform of comparison and the match is one by definition. 



In order to have a qualitative understanding of why the match behaves as in Fig. [H 
we study the p per mode for the A-series. This is only a qualitative estimate because the 
relative fraction of modes present in the signal will depend on the relative spin-weighted 
spherical harmonics values at the particular angle. In practice the error induced by 
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ignoring the mixed terms that are important in constructing p{6, 0) from the modes is 
less than 20%, as the relative overlaps of the significant modes from both the A and 
B-series of data are of this order. Since the p of the i = m = 2 mode is much larger 
than the p of the other modes, we plot the ratio, p(£, m)/p(2, 2) in Fig. [2l The upper 
left plot corresponds to a system of mass SOM©, the upper right to lOOM©, lower left to 
2OOM0 and lower right to SOOM©. Across all the masses sampled, the ratio of the p for 
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Figure 2. A series of plots are slrown witli the ratio of p per mode to the p of the 
£ = m = 2 mode versus the initial spin of the black holes. This is done for series-A. 
Each plot refers to the calculation for a different total mass of the binary. Starting 
from the upper left and moving right and then down, we have 50AIq, IOOMq, 200Mq, 
and 3OOM0 on the lower right. 



each mode grows with increasing spin. This is especially true for the m = 1 and m = 3 
modes which are suppressed at low a. The m = 2 and m = 4 increase slightly with a. 
While the £ = 2, m = 1 mode is the next mode dominant mode after the i = m = 2 
mode for the high-spin regime at low masses, the i = m = 4 mode is secondary for the 
entire a range at higher masses. For low spins, the waveform is entirely dominated by 
the m = 2 modes. The linear-like growth of the odd-m modes with spin magnitude is 
expected from PN expressions like Eq(l)-(4) in [2T] . 

B-series: The minimax match versus the initial angle for the B-series is presented in 
Fig. [3l Each line represents a choice of total mass, with SOMq the top most line and 
3OOM0 the bottommost. The match was computed by setting one waveform to ^ = 
and the other to ^ = 7r/2.4. We find that the variation of the match across initial 
angle for the given spin of a = 0.6 does not change more than about 2%. The variation 



Spinning Binary Black Holes 



6 



1 





^bm^ 


0.98 




100M. 


0.96 








0.94 




0.92 


250M^ 


0.9 


300M^ 


0.88 





0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 

t> 

Figure 3. The minimax match versus initial angle d for the B-series. Each curve 
represents a particular choice of total mass, SOM© at the top with each successively 
lower line a higher mass. To compute the match, we used one waveform with 6 ~ Q 
and the other at 9 = 7r/2.4 radians for a given a and §. 

amongst different total masses is more dramatic, dropping down below a match of 0.9 
for most in the angles at a mass of SOOM©. At that large mass range, the ringdown is 
contributing significantly to the signal, and differences in the modes, like the £ = m = 4 
mode, begin to make important contributions. These BBH configurations settle down 
to a final black hole with a spin of a = 0.62. 

In Fig. m we once again investigate a qualitative interpretation of the match through 
the p as plotted versus the initial angle, for each mode. The upper left plot corresponds 
to a system of mass 5OM0, the upper right to lOOM©, lower left to 2OOM0 and lower 
right to SOOMq. Across the mass scales sampled, as d increases, the signal in odd-m 
modes decrease. In the non-precessing case, the strength of the odd-modes is expected 
to vary with the z-component of the spin [21]. Since in this series of runs, the spins 
precess about the z-axis and d remains nearly constant, the relative strength of the 
modes show similar trends as the non-precessing case. At higher masses, where the 
ringdown dominates the signal, the £ = m = 4 mode contributes a large portion of the p 
of the total signal, over 20% for a BBH of SOOM©. This is in part due to the enhancement 
of the ringdown signal which occurs when the frequency of the higher modes lies around 
the detector's sweet spot. It is interesting to note that in the "superkick" the spread 
of the modes is reduced in p compared with the parallel configuration aX d = Q. The 
decrease of the odd m modes is expected from the PN expressions. For example, the 
£ = 2, m = 1 will be suppressed when the spins lie in the orbital plane for equal-mass 
black holes as discussed in [21] . 

These waveforms are the solution to the BBH coalescence as expressed by general 
relativity with errors arising from several sources, see [24j. In paper I, we analyzed the 
effects of resolution on the matched filtering technique and found that for the resolutions 
used to compute the waveforms studied in the A-series, the largest error would be ±0.02 
in the match, although that is only for the a = 0.8 case, and is typically smaller. The 
waveforms in the B-series have comparable errors, i.e. the resolutions, wave extraction 
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Figure 4. A set of plots is sliown with the ratio of p per mode to the p of the 
£ = m ^ 2 mode versus the initial angle for series-B. Each plot is the ratio computed 
for a different total mass of the binary. Starting from the upper left and moving right 
and then lower left and right, we have SOM©, lOOMg, 2OOM0, and SOOM© on the 
lower right. 



and other numerical techniques were the same in computing both series of waveforms 
as discussed in P [151 IB]- For reference, the typical resolution on the finest grids were 
h = M/35.2 where M is the total mass. 

4. Discussion and Conclusion 

In this paper, we investigated the contribution of individual modes to p from the last 
orbits, merger and ringdown of an equal-mass, spinning binary black coalescence. In the 
A-series, the spins were kept parallel/anti-parallel to the direction of the orbital angular 
momentum, but the magnitude of the spins varied. In the B-series, the magnitude was 
kept fixed to a = 0.6, but the initial angle the spins make with the orbital angular 
momentum varied. 

In paper I, we investigated the match between a waveform from the A-series 
containing only the i = m = 2 mode and a waveform of a sum of modes. There we 
found strong dependence on the match with spin, with the i = m = 2 waveform failing 
to match to spinning waveforms especially for spins equal to and greater than a = 0.6. 
We did a similar study here for the B-series, comparing two waveforms of a = 0.6 at 
various {}. We found, despite the variation of the p versus {}, the match had a much 
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greater dependence on mass than the initial angle. The inclusion of modes was much 
more important to templates of higher mass, where the merger and ringdown dominate 
the signal, than at lower masses. This importance will be more evident in matches using 
unequal-mass and spinning waveforms with larger spin as well as waveforms with more 
cycles. 

To qualitatively understand the matches, we conducted a multipolar analysis of 
the modes in each waveform and calculated the ratio of the p of each mode versus 
C. = m = 2. For the A-series, the p per mode increased as the magnitude of the spins 
increased at every total mass as seen in Fig. [2l The odd-m modes increased from almost 
no contribution at low spins to a 10% contribution at larger spins. At a given spin, the 
£ = 2, m = 1 mode dominated the p at low mass, but the £ = m = 4 mode's ratio 
to ^ = m = 2 grew with increasing mass. For the B-series, Fig. HI we found that the 
diversity of contributing modes decreases with increasing angle, except for the £ = m = 4 
and £ = 3, m = 2 modes which remain relatively constant across d for a given mass. 
As in the variation with a, at low total binary mass, the secondary signal is the ^ = 2, 
m = 1 mode, but at higher masses the £ = m = 4 and the I = m = 2 modes are 
stronger. As anticipated, the £ = 2, m = 1 mode decreases to zero as the initial angle 
moves to lie parallel to the orbital plane. For both the series of runs, the variation of 
the signal in different modes is consistent with the expectation from PN [21] . 
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